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Compact stars consisting of fermions with arbitrary masses and interaction strengths are studied 
by solving the structure equation of general relativity, the Tolman-Oppenheimer-Volkoff equations. 
Scaling solutions are derived for a free and an interacting Fermi gas and tested by numerical calcu- 
lations. We demonstrate that there is a unique mass-radius relation for compact stars made of free 
fermions which is independent of the fermion mass. For sufficiently strong interactions, the maxi- 
mum stable mass of compact stars and its radius are controlled by the parameter of the interaction, 
both increasing linearly with the interaction strength. The mass-radius relation for compact stars 
made of strongly interacting fermions shows that the radius remains approximately constant for a 
wide range of compact star masses. 
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I. INTRODUCTION 



Compact stars, white dwarfs and neutron stars, are one of the most enigmatic astrophysical objects (for an in- 
troduction to the physics of compact stars we refer to the excellent textbooks Q, |^ and to recent pedagogical 
papers 0,11, IE SIB])- The first successful description of compact stars was proposed by Fowler in 1926 Q, who first 
realized that Fermi-Dirac statistics is responsible for the high degeneracy pressure which holds up the star against 
gravitational collapse. Shortly afterwards Chandrasekhar applied effects of special relativity to the equation of state 
(EoS) of a degenerate Fermi gas and demonstrated the existence of a maximum mass for such stars beyond which 
they are unstable against collapse: the famous Chandrasekhar mass limit, Mch ~ IAMq p^ . 

In 1932 Chadwick discovered the neutron and Heisenberg presented his theory of isospin for nucleons suggesting that 
the neutron is a fermion with spin-i like the electron. The ideas of Fowler and Chandrasekhar were then immediately 
extended to the case of degenerate neutrons and a new form of compact stars, neutron stars, was predicted by Landau 
|llj |. The first numerical calculations for a neutron star within general relativity were performed by Oppenheimer 
and Volkoff in 1939 They computed a maximum stable mass of O.TSMq for a free Fermi gas of neutrons beyond 
which the star would be unstable and collapse into a black hole. The maximum mass limit for neutron stars is now 
well known as the Oppenheimer- Volkoff mass limit. Contrary to the case of the electron gas in white dwarfs, there 
are sizable contributions to the mass limit for neutron stars due to strong interactions between the neutrons. These 
interactions can be well modelled by an effective repulsive potential which can increase the maximum stable mass of a 
neutron star to about 3Mq . The presence of hy pero ns in the core of neutron stars leads to a softening of the equation 
of state and a reduction of the maximum mass |l3, [ij, [T^ [1^, Il3 ■ 

As of today, new massive fermions are known within the standard model and many more have been predicted, in 
particular also as candidates for dark matter. In the year 1964, Gell-Mann and Zweig proposed the idea of quarks, 
putting forward the notion that neutrons and protons are composed of quarks. Ivanenko and Kurdgelaidze [l^ studied 
a quark core in massive compact stars and Itoh |19| calculated the mass-radius relation of a quark star. If a compact 
star consists of quarks only, including besides the light quarks also the strange quark, they are dubbed strange stars 
[20ll2lll2^ . Those exotic objects might be bound by strong interactions only, contrary to ordinary neutron stars and 
white dwarfs which are bound by gravity. The physics of compact stars with a quark core and strange stars is now 
a field of active research (for a recent review we refer to [23]) and has found its place in modern textbooks |3i0- 
Besides quarks, other fermions in the form of massive neutrinos are now well known to exist in nature. New types 
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of fermions predicted in extensions of the standard model include the supersymmetric particles, the neutralino, the 
gravitino, and the axino, which are also candidates for dark matter (for a review see e.g. p^*). 

Now the idea put forward by Fowler many years ago could be used to speculate on compact stars made out of exotic 
fermions, stabilised by the degeneracy pressure in the same way as in the case of neutron stars and white dwarfs. The 
present investigation is dealing with the following questions: What is the maximum stable mass of compact stars as 
a function of the fermion mass at zero temperature? What will happen if a repulsive interaction is included in the 
equation of this fermionic matter? 

The paper is organised as follows: in section ^ we briefly recapitulate the structure equation for compact stars, 
the Tolman-Oppenheimer-Volkoff (TOV) equation. In section HTll we discuss general scaling solutions for compact 
stars with an arbitrary equation of state, in particular for a free Fermi gas and for an interacting Fermi gas. We show 
that Landau's estimate for the maximum mass and the corresponding radius is an exact scaling solution of the TOV 
equation for a free Fermi gas. Section IWI is devoted to our numerical solution of the TOV equations. The equation 
of state for interacting fermions as well as extended scaling solutions of the TOV equations are discussed. We show 
that there is one unique mass-radius relation for compact stars made of free fermions if it is properly rescaled by the 
Landau mass and the Landau radius. We also find that for strongly interacting fermions the maximum mass and its 
radius are controlled by the interaction terms. The mass-radius relation changes such that there is a constant radius 
for a wide range of masses. Finally, in section^we summarise and discuss our findings. 

II. THE STRUCTURE EQUATION FOR COMPACT STARS 

Throughout this paper we will be dealing with models of compact stars where effects of general relativity are 
included for a consistent description as in the case of ordinary neutron stars. The typical mass of a neutron star is of 
the order of ~ IM© with typical radii of about 10 km, i.e. lO~^i?0. Hence, the gravitational potential on the surface 
of neutron star will be 10^ times stronger compared to that of the Sun. Moreover, the corresponding Schwarzschild 
radius = IGMjf? is about 3 km in this case. Under such conditions, the curvature of space-time can not be ignored 
and general relativity is needed to describe the structure of such compact objects. 

In order to find the structure of space-time created by the presence of a compact star one needs to solve the Einstein's 
field equations. For simplicity we assume that the metric is spherically symmetric and static, i.e. the Schwarzschild 
metric. In addition, the energy momentum tensor is assumed to be that of an ideal fluid, 

T^i^ ^ pg,,,^ + {p + p)u^,u^ , (1) 

where p and p denote the pressure and energy density, and C/^ the four velocity of the fluid. Using the Einstein's field 
equation and the condition for hydrostatic equilibrium, Uf^ = (liO)j one arrives at the following equations describing 
the structure of a compact star 




which are just the Tolman-Oppenheimcr-Volkoff (TOV) equations 0, (note, that throughout the paper, we 

are using natural units by setting h — c — 1). The detailed derivation of the equation (0) can be found in standard 
textbooks (see e.g. [H S 123 ) ■ This is the same equation of hydrostatic equilibrium as in the case of Newtonian 
gravity just modified by three correction factors (in the r.h.s.) due to effects of general relativity. The equation Q 
simply defines the quantity M{r), the amount of energy contained within the radius r. 

The unknown functions in eq. ([SJ and eq. Q are p{r), p{r) and M{r). For a given equation of state, relating p and 
p, appropriate initial and boundary conditions are needed to solve the above set of equations. The radius of the star, 
R, is found by using the condition that the pressure vanishes at the surface of star. The mass Af(0) must be zero at 
r = and M{R) gives the total mass of the star at r = i?. The central pressure is calculated from the equation of 
state once the central energy density p(0) — po is given as the initial condition. 

III. SCALING THE TOV EQUATION 

It is easy to see, that the TOV equation contains two essential dimensional quantities, G which can be conveniently 
expressed in terms of the Planck mass as G = M~'^, and the fermion mass m/ which characterises the equation of 
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state. In this section, we show how the TOV equation can be transformed to scale-independent variables composed 
of Alp and ruf. There are several reasons for such an approach. One reason is that the computational treatment of 
differential equations benefits from a dimensionless format. The other reason is that a scaled equation needs to be 
solved only once. As soon as the general solution is found one can just rescale it by appropriate (dimensionful) factors 
to get the result for specific (astro-)physical cases. The TOV equation can be scaled in the following ways. 



A. Landau's argument for deriving the maximum mass of compact stars 

Landau presented a very elegant argument for deriving the maximum mass of a compact star llj (for details see 
also the treatises in [2,, 29, ,3QJ). He used only Newtonian gravity, special relativity and Fermi-Dirac statistics to 
estimate the maximum mass and the corresponding radius of a compact star. 

For a star made of free fermions at zero temperature the Fermi momentum fc^^^ is related to number density n by 
the following relation (below we suppress the dependence on the statistical degeneracy factor assuming g = 2) 

_ _fc|^ _ N 

" 3^2 - 47r/3i?3 ' 

where N is the total number of fermions in a star. For simplicity, here we consider a star of uniform number density 
and of radius R. Solving for kp, one gets 

-^(?)"^- 

The total mass of the star is given solely by the vacuum fermion mass, to/, as 

M = rufN . (6) 
Now, let us consider a fermion on the surface of the star. Its energy is given by 

where the first term gives the gravitational energy and the second term comes from the kinetic energy of an ultrarel- 
ativistic fermion on the top of the Fermi distribution. For small values of R and negative energy E, the gravitational 
attraction overcomes the degeneracy pressure causing a collapse. For positive E the degeneracy pressure exceeds the 
gravitational attraction and the star will expand until the particle density drops so much that kp ~ to/. For the 
non-relativistic gas the kinetic energy per particle is 3fc|,/(10TO/), i.e. it changes with the radius as R~^. This means 
that gravitation will finally prevail and the expansion will be terminated. As a result, a stable minimum will develop 
in E{R). 

A good estimate of the maximum possible number of fermions in a compact star can be obtained by considering 
the limiting case E — 0, when the gravitational energy is exactly equal to the Fermi energy of the degenerate Fermi 
gas. Then we can express the maximum number of fermions as 



97r\ M3 



4 J my 

The maximum mass of the star is obtained from the relation M = irifN, hence 

p 

'a 



3 • (8) 



my 

An estimate for the corresponding radius of the maximum mass star can be obtained by assuming that the kinetic 
energy of the fermion on the surface is equal to its mass, i.e. at the border of becoming relativistic, fc^? « to/, which 
gives for the minimum radius of the star 
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For a neutron star, with the fermion mass taken to be that of the neutron m„ « 1 GeV, the above relations give 
Mmax ~ 1.63Mq and Rmin ~ 2.41 km which is a reasonable estimate. Note, that Landau's argument can be well 
applied also for white dwarfs. Here, one has to take care of the fact that the mass of the white dwarf is determined 
by the nucleon mass while the degeneracy pressure is provided by the electrons. Therefore, the maximum mass for 
white dwarfs, the Chandrasckhar mass limit, turns out to be similar to the one for neutron stars, while the radius 
increases by the ratio of the nucleon and the electron masses, i.e. by about a factor 2000, to about 4000 km (all these 
values have to be corrected for the charge to mass ratio of nuclei, see e.g. [30j). 



B. The equation of state in dimensionless form 



The equation of state for a free gas of fermions at zero temperature p{p) can be calculated via explicit expressions 
for the energy density and pressure: 



'-2 /mj 



k^dk 



1 
3^ 



247r^ 



{2z^ + z)(l + z^) 2 - sinh^^z) 

dk 



_ 4 / 



(11) 



fc2 



(2z3-3z)(l + z2)2 +3sinh"i(z) 



4 / 

rajp 



(12) 



using natural units of and defining the relativity parameter z = kp/irif. 

We introduce now the following dimensionless quantities for the mass and the radius of the star. 



M Ml 
= — with Ml = ^ 
Ml mj 



and r — with Rl = — ir 

Rl mi 



(13) 



where Ml and Rl denote the maximum mass and corresponding radius as given by Landau's arguments (note that 
Rl is equal to half the Schwarzschild radius). Using the above definitions and substituting them into eq. Q and 
eq. Q with G = Mp'^, one obtains the following dimensionless form of the TOV equations: 



d? 
dM' 
~d?~ 



M'p' 



M' 



2M' 



(14) 
(15) 



In general, the equation of state can not be expressed in a simple polytropic form p' cx p'^ with a constant 7. However, 
this can be done for two limits of the relativity parameter z. In the non-relativistic limit, z ^ 1, we get a polytropic 
law with 7 = 5/3, i.e. 



/ /5/3 
p (X p 



(z«l) . 

In the ultra-relativistic limit, z ^ 1, the equation of state approaches a polytrope of 7 1, i.e. 



(16) 



(17) 



As the equation of state is a function of the relativity parameter z, one can compute it parametrically in a tabular 
form for a desirable interval of z. Fig. ^ depicts the resulting dimensionless pressure versus the dimensionless energy 
density in a double logarithmic plot. The curve exhibits two different slopes for small and large values of p' . The larger 
slope at small values of p' shows the non-relativistic regime, the smaller slope at larger values of p' the relativistic 
regime as expected. 

For the case of white dwarfs, the following changes have to be made to the equation of state, eq. Hll|l and eq. I|12(l : 

1. As the pressure is given by the degeneracy pressure of electrons, the expression for the pressure remains the 
same where m is the mass of the electron. 
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FIG. 1: The equation of state, the dimensionless pressure p' versus the dimensionless energy density p' for a free Fermi gas at 
zero temperature. The two limits, the non-relativistic one for z and the ultra-relativistic one for z ^ 1, can be clearly 
seen. 



2. The energy density is a sum of two terms, one coming from the mass density of protons and neutrons while the 
other from the kinetic energy of electrons. 

For a white dwarf consisting of nuclei with the atomic mass A and atomic number Z the mass density of protons and 
neutrons reads (small effects from the nuclear binding energy are disregarded here) 

A , , 

pi=n- niN ■ — (18) 

where due to charge neutrality the number density of protons n is equal to the number density of electrons which is 
given by eq. (0J). The total energy density is now written as 



P = Pi 



-i- / xWl + x^dx . (19) 

7^ Jo 



The second term here gives the relativistic energy density of electrons (note that here z — kp/me). The resulting 
equation of state can be scaled in the same way as done before. Defining the dimensionless pressure and dimensionless 
energy density as p' — p/m^ and p' = p/m^, we can rewrite the above equation in the dimensionless form 



, TUN A 1 2 n ^1 

6-1^^ me Z TT^ Jq 



(20) 



The pressure p' is given by eq. (|12|l where mj = rrie- In the non-relativistic case, z <C 1, the equation of state simplifies 
to eq. (|16|l . On the other hand, in the relativistic limit 1 ^ z ^ 6 • 10^, when the first term on the r.h.s. dominates, 
the equation of state for white dwarfs becomes 



p'cx//' . (21) 

The proportionality constant in both limits, non-relativistic and relativistic ones, depends on A and Z. Note, that the 
special case for white dwarfs emerges due to the presence of oppositely charged fermions which are neutralising each 
other. The heavier one determines the energy density, while the lighter one the pressure. For a pure charge-neutral 
Fermi gas, as we will consider later, the ultra-relativistic limit will be realized, as for the case of neutron stars. 
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C. Scaling solution for self-bound compact stars 

Another scahng solution is well known [l^,!^ for the special form of the equation of state (see also 0| and references 
therein) : 

p = s{p-pa) (22) 

where s and po are constants. Note, that this equation of state has a vanishing pressure at a finite energy density 
p — Pq. This property allows for the existence of self-bound balls of any size, not necessarily of astronomical scale. 
These balls are stabilised by other interactions not by gravity as in the case of neutron stars. 

The scaling relations p' = p/po, p' = p/po, r' — ^/po{r/Mp) and M' = yjrp^[M/Mp) change the TOV equation 
to a dimensionless form. One can then solve this equation numerically. If the mass and the radius are known for 
some particular value of po, then for some other value p'^ the radius and mass will be R{p'q) = \J p'^/ poR{po) and 
^{p'o) = V p'o/PoM{po), respectively. Hence, both the mass and the radius scale with Xj ^fp^. 

The above equation of state is actually the one of the MIT bag model often used for describing cold and massless 
(strange) quark matter. The corresponding self-bound compact stars are dubbed strange stars |23,0,|23. Most 
studies of quark stars and strange stars utilise the equation of state derived from the MIT bag model. Ignoring effects 
from a finite quark mass, this model gives the equation of state 

p=i(p-4B) , (23) 

where B is the bag constant. Interestingly, the equation of state for an interacting cold gas of massless quarks within 
perturbative quantum chromodynamics can be approximated by the same form of the equation of state ^31j . Standard 
values for the MIT bag constant are around B^l^ — 145 MeV as follows from fits to hadron masses, which results 
in maximum masses of about 2S)Mq at a radius of about 11 km pol l2ll Is^ . which are actually very close to the 
ones of realistic neutron star models. It is worth noting, that these values of B are obtained in fits including the 
chromomagnetic interaction. In simplified versions of this model, disregarding the interaction effects, higher values, 
B^l^ ~ 200 MeV, are needed to preserve the stability of normal nuclear matter |33| . 



D. Scaling in a general case 

Consider the TOV equation for the pressure, eq. We make the observation that the three correction factors 
from general relativity are already in a dimensionless form. Therefore, from the first factor (1 -\-p/p) one can scale 
pressure and energy density by a common factor eo as p' = p/co and p' ~ p/fo, respectively. Similarly, one defines 
a dimensionless mass and radial coordinate via M' = M/a and r' = r/b. Plugging these definitions into the second 
factor (1 + ^irr^p/M) one arrives at a dimensionless number b^eo/a which we equate to one. Similar reasoning for the 
third factor (1 — 2GM/r)~^ results in a/{Mpb) as the dimensionless number which again is equated to one. So we 
get finally the following scaling conditions 



a M^b 
Solving for a and b from the above equations one gets the following expressions 



1 and T72l = l • (24) 



and 5 = — £: . (25) 



For a free gas of fermions at zero temperature, we know from the direct calculation that eo = (see eq. Illll . 
Substituting this value of eo into the above equations one finds a = Mp/rnj: and b = Mp/rnj: which are exactly the 
same scaling factors as used originally by Landau. In the case of the MIT bag equation of state, the scaling factors 
are given by eqs. with eo = po = 4i3. 

We note in passing that according to Landau's argument it is sufficient to incorporate special relativity and New- 
tonian gravity to obtain the maximum mass Mmax and the corresponding minimum radius Rmin of a cold compact 
star made of fermions. Interestingly, the same dimensional forms for the maximum mass and minimum radius are 
found when using dimensional reasoning applied to the full TOV equation of general relativity. 
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FIG. 2: Mass-radius relation in Landau units obtained by solving numerically the dimensionless TOV equation for a range of 
dimensionless central densities from p' — 10~* to p' — 10. Note, that the curve does not depend on the mass of the fermions 
forming the compact star. 



IV. NUMERICAL SOLUTION OF THE TOV EQUATION 
A. Stars made of free fermions 

For a free gas of fermions at zero temperature, the equations to be solved are the dimensionless TOV eq. (|14ll and 
eq. H15(l along with the dimensionless equation of state given in parametric form by equations and (|12|l . One 
starts vifith the dimensionless central pressure or energy density, then solves the TOV equations from the center of 
the star to the surface where the pressure becomes zero. The corresponding radial distance defines the radius of the 
compact star. 

Since the dimensionless equation of state is given as a function of parameter z, both quantities p' and p' are 
calculated for various values of z and expressed in a tabular form. To calculate p' for a given p' and vice-versa, 
we use a simple linear interpolation. The dimensionless TOV equation is solved using a fourth-order Runge-Kutta 
algorithm. To plot the M' versus B! curve, M' and B! are calculated for various dimensionless central densities. 
The step size of the dimensionless radius lies between 0.01 (for very low dimensionless central density) and 0.001 (for 
very high dimensionless central density) to have at least 2000 points for one star configuration. The step size in the 
dimensionless central density is adjusted to have 500 points between lO^'^ and 10. The final dimensionless mass-radius 
relation is plotted in Fig. El 

From the curve in Fig. [21 one notes that for a small dimensionless mass the dimensionless radius is large. This 
behaviour occurs for very small dimensionless central densities. The gravitational attraction is small in this case 
making the dimensionless radius large. As the dimensionless central density increases, the mass of the star increases, 
too. This leads to a stronger inward gravitational pull and smaller dimensionless radii. Thus, as we increase the 
dimensionless central density the mass increases while the radius decreases. A maximum mass is reached for the 
dimensionless mass M^^^ = 0.384 at the dimensionless radius i?^j„ = 3.367. The presence of the maximum in the 
curve is generic and expected for an arbitrary fermion mass. The reason is that the energy density which generates 
the gravitational pull inwards has to be balanced by the outward fermionic pressure. However, the rate of change 
of the pressure with energy density is related to the speed of sound which is bounded by the speed of light. This 
speed limit puts a bound on the pressure increment with respect to changes in energy density. Thus an increase 
of the central energy density results in an increased gravitational attraction which cannot be compensated by the 
corresponding additional pressure, that leads eventually to a maximum mass limit ^3^3- easy to see Q, that 
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FIG. 3: The above graph shows the M' versus R' curve on a logarithmic scale for large radii. 



star configurations on the right hand side of the maximum in Fig. [5] are stable whereas those on the left hand side 
are unstable. The maximum in the curve proves the fact that there exists a maximum stable mass for a fermion star, 
beyond which the stars are unstable and collapse. 

Fig. O shows M' versus R' on a double-log scale. For large values of the dimensionless radius i?', the mass M' and 
the radius R' follow the relation 

M' ■ R''^ = constant (26) 

A linear fit to the curve for large R' gives (3 = 3.0001 « 3 and the constant is 90.97 w 91. Putting these values in the 
above relation we get 

M' ■ R'^ « 91 (27) 

The above relation can be also derived analytically by considering a non-relativistic degenerate gas of free fermions 
and Newtonian gravity. Compact star configurations with a polytrope form of the equation of state, p ~ p''', are given 
by the Lane-Emden function [23, 0| . A non-relativistic gas of fermions has a polytrope of 7 = 5/3 and obeys the 
Lame-Emden equation of index (7—1)^^ = 3/2. For a general polytrope, the relation between M', R' and the central 
density p'(0) are given by 

M' oc p'(0)(3^-4)/2 and i?' cx ^'(O)^^-^)/^ =^ M' oc i?'^'^-"^/^^-') (28) 

and for a polytrope of 7 = 5/3 one finds M' oc i?' ^ , exactly as in eq. H27|l . 

Now we can clearly see the advantages of solving the dimensionless TOV equation. To get full dimensional results 
for a given fermion mass one should simply rescale the dimensionless mass and radius with factors a = Mp/rrvj and 
b ~ Mp/nij, respectively. For the fermion mass mf = 1 GeV, these factors are a — 1.632Mq and b = 2.410 km. 
Therefore, for a point on the curve in Fig. |2] with the dimensionless mass M' and the dimensionless radius R', the 
actual mass and radius for a fermion star will be given by 

M = 1.632 Mq-M' ■ (^-^^\ (29) 
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Fermion mass 


Mmax(M0) 




comment 


100 GeV 


10-* 


1 m 


neutralino star (cold dark matter) 


1 GeV 


1 


10 km 


neutron star 


1 GeV/0.5 MeV 


1 


10^ km 


white dwarf 


10 keV 


10'° 


10" km 


sterile neutrino star 


1 keV 


10'2 


10" km 


axino star (warm dark matter) 


1 eV 


10'« 


10'^ km 


neutrino star 


10-2 eV 


10^2 


lO^^* km 


gravitino star 



TABLE L Maximum masses Mmax and radii Rmin for various cold compact stars made of a free Fermi gas 



and 



R = 2.410 km -R' ■[ . (30) 

ruf J 



Accordingly, the maximum mass and minimum radius of a fermion star are obtained for M' = 0.384 and R'min — 3.367 



Mma. = 0.627 Mq ■ ( (31) 



TO/ 



and 



i?™„ = 8.115 km • (l^) (32) 

One can use the above relations to calculate the masses and radii of different fermion stars. For example, a neutron 
star with to/ = to„ — 939.6 MeV in eq. (|31|l and eq. H32|l . has the maximum mass Mmax = 0.710 Mq at a radius of 
Rmin — 9.192 km, which match very well the original results of Oppenheimer and Volkoff [l^ . 

Next, we consider compact stars built of other fermions utilising the curve shown in Fig. El Besides the nucleon and 
electron, fermions as neutrinos and quarks are well established in the standard model. Moreover, other fermions such 
as neutralinos, axinos and gravitinos are predicted in supersymmetric extensions of the standard model. For ordinary 
neutrinos, the mass has been recently constrained to the r ang e from 50 meV to about 1 eV from the measurements of 
neutrino oscillations and cosmological parameters (see [35ll36j |). Hypothetical sterile neutrinos can have typical masses 
in the keV range (see e.g. [s^). For the supersymmetric particles, neutralinos, gravitinos and axinos, there are large 
uncertainties concerning their mass ranges. The most likely mass for neutralinos is usually considered to be around 
100 GeV but lighter neutralinos with a mass of < 50 GeV are discussed in the literature, too [s^- Gravitinos can be 
heavy or very light in some models, the latter mass range can extend from 10^^ eV to 1 keV, bounded by constraints 
from big bang nucleosynthesis and the critical density (see e.g. the discussion in |39j|). Axinos, the supersymmetric 
partner of the axion, were introduced first as a possible warm dark matter candidate with a mass in the keV range 
pol |. but they are discussed now also as a cold dark matter candidate with much higher masses Compact stars 
could be also formed from composite dark matter candidates, similar to white dwarfs, e.g. with the heavy charged 
stable leptons proposed in 42]. Note, that the preceding discussion is just to motivate interesting mass ranges for 
fermions to be considered in the following. There are of course constraints for dark matter candidates per se which 
are more severe than the ones mentioned above, see e.g. the recent review on dark matter candidates in p3 |. 

In principle, cold compact stars can be formed out of these exotic fermions, too. Actually, compact stars made of 
massive neutrinos have been introduced by Markov and calculated within general relativity by Gao and Ruffini 
|44|. Using the above mentioned masses for these fermions, we have calculated their typical (i.e. maximum) masses and 
corresponding radii which are listed in Tabled A compact star made of non-interacting neutralinos has a maximum 
mass of about 10-* solar masses with a radius of about one meter. Warm dark matter, fermions with a mass in the 
keV range, can form compact objects with galactic masses and a radius of about one light year. In ref. [islli^ sterile 
neutrinos with a mass of m/ ~ 50 keV were proposed to explain the dark massive object at the center of our Galaxy, 
as an alternative to a supermassive black hole. Compact stars made of dark matter interacting with a scalar field 
have been also considered in [47„ .48] . Interestingly, if we take the gravitino mass to be of the order of 10-^ eV, then 
the corresponding gravitino star has the mass and the radius of our universe. Note, that the above masses and radii 
are similar to the free-streaming mass and radius scales which are of importance for large-scale structure formation. 
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B. Stars made of interacting fermions 

1. Equation of State for Interacting Fermions 

The previous section was devoted to an idealised case of non-interacting fermions. In a more realistic consideration, 
the inter-particle interactions must be included, too. In the following, we address the possible impact of interactions 
on the global properties of fermion stars. 

Consider the simplest model of two-body interactions between the fermions. In a lowest order approximation the 
interaction energy density is proportional to n^, where n is the number density of fermions. To have the correct 
dimensionality this term can be written as pint — n'^/m'j where mj represents the energy scale of the interaction. The 
corresponding contribution to the pressure is: 



P 

oV 



lid Pint /n) 



N,T=0 



dn 



m. 



(33) 



Therefore, in this approximation the energy density and pressure acquire an additional term /mj. The interaction 
must be repulsive so that an increase in the number density increases the pressure and energy density. The scale m/ 
can be also interpreted as the vacuum expectation value of the Higgs field of the interaction. For weak interactions, 
the interaction strength is just given by Fermi's constant or the vacuum expectation of the Higgs field v generating 
the masses of the W and Z bosons, Gj,/V2 = l/{2v^) (293 GeV)-^, so that m/ ~ 300 GeV. Correspondingly, the 
strength of low energy strong interactions, quantum chromodynamics (QCD), is controlled by the pion decay constant 
1//^ in chiral perturbation theory, with = 92.4 MeV being the vacuum expectation value of the sigma field (for 
an introduction to chiral symmetry see e.g. ^3)- quantum hadrodynamics, the expressions for the energy density 
and pressure are exactly as given above for a repulsive interaction mediated by a vector meson with an interaction 
strength of 5f,Ar/ which is quite close to 1//^ for g^jy — 13 and — 780 MeV (see e.g. H^l). Hence, for strong 
interactions the typical interaction mass scale is mj ^ 100 MeV. In a more realistic approach, an attractive scalar 
interaction should be included in addition to a repulsive vector interaction, too. This kind of approach is widely used 
for baryonic matter 50] as well as for quark stars (see e.g. HJ). In dimensionless variables the energy density and 
pressure can be written as 
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(34) 



(35) 



where z is the dimensionless Fermi momentum and y = mf / mi the interaction strength. 

What are the values to be taken for the interaction strength yl For a realistic neutron star, neutrons interact 
strongly with mi ~ 100 MeV, as outlined above. With a neutron mass oi mf ^ 1 GeV one arrives at y '--^ 10. 
Neutrinos interact weakly and with a neutrino mass of about 1 eV (for sterile neutrinos m/ ^ 1 keV), one finds 
y ^ 10~^^ (for sterile neutrinos y ^ 10~^). For neutralinos with a mass of 100 GeV and weak interactions, y ~ 1/3, 
for strongly interacting neutralinos y ~ 10'^. We find that for small values of y the mass-radius relation remains 
almost unchanged because the equation of state is dominated by the kinetic terms. Any change in the equation of 
state can only occur for ?/ > 1, i.e. when interaction term starts dominating the equation of state before the fermions 
are becoming relativistic. Hence, it is sufficient to take y to be in the range 10"^ to 10"^. Below y = 10"^ one will 
hardly observe any change in the mass-radius relation. 

From the equations for p and p one finds that the common factor mi can be taken out, so that general scaling 



arguments result in eo = m^^. Therefore, the dimensionless forms oi p and p are p' = p/m'^f and p' 



M/a and R' = R' /h with M^/ 



f 



p/m'j: and 



mj and b = Mp/mj, exactly 



the corresponding mass M and radius R are M' 
Landau's mass and radius as before. 

Note, that for each value of y there is a different equation of state and two different regimes exist: z ^ 1, the 
non-rclativistic limit, and z ^ 1, the relativistic limit. For small y <C 1, the equation of state will be that of an 
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ideal Fermi gas and for large y 3> 1 the equation of state will be mostly determined by the interaction term, unless z 
becomes small enough so that the ideal gas term becomes the dominant ones. 

The dimensionless TOV equations (|14|) and (|15|l are solved by making use of the dimensionless equation of state, 
eqs. (|34ll and (|35|) . The equation of state depends on two parameters, z and y. To solve the dimensionless TOV 
equations for particular values of y, we first construct the dimensionless equation of state in a tabular form for 
different values of z. Linear interpolation is used to find p' corresponding to a particular dimensionless density p' 
and vice-versa. As before, the dimensionless TOV eq. H14|l and eq. H15|l are numerically solved using a fourth order 
Runge-Kutta algorithm. To plot the graph of M' versus B! for various y, we take 100 equally spaced values of y lying 
between 10~^ and lO'^. For each value of y, the M' versus E! curve is computed for 50 equally spaced dimensionless 
central densities. The program adjusts the step size in the dimensionless radius to have about 2000 points for each 
star configuration. Fig. 0] shows the final logp' versus logp' curve for a range of values of z from 10^^ to 100 for six 
different values of y. 

From Fig. 0] one makes the following observations: 

1. For large values of p', the different equations of state merge to one line, except for the curve of the cases y = 
and y = 0.01 which stay slightly below the other ones. 

2. The transition point from the relativistic curve to a steeper one shifts to lower densities for increasing interaction 
strength y. 

3. For very small values of p' , the slope of the curve is approaching the slope for a non-interacting ideal gas. The 
higher the interaction strength y, the lower p' must be taken to reach the non-interacting limit. 

These observations can be explained in the following manner. The equations H34|l and H35|l are parametric in z and 
y. In the non-relativistic limit, z <C 1, they are reduced to 

and 
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In the relativistic limit, z ^ 1, one finds 

+ ^^^^ 

and 

If the interaction strength is small, y <C 1, the interaction terms in eqs. H36(l and can be ignored for z ^ 1 and 
one recovers for small densities the non-interacting case. For increasing y, the interaction terms become more and 
more important. Let us examine the ratio of the interaction term to the free gas term in more detail for the different 
cases. We define the two ratios Tp{z ^ 1) and Tp{z <^ 1) for the dimensionless energy density and pressure eqs. H36|) 
and IpTf)! as 

Tp(z < 1) - y^z^ and rp(z < 1) - y^z (40) 

The corresponding ratios Tp{z ^ 1) and Tp{z ^ 1) for the case z ^ 1 follow from eqs. H38|l and (|39(l . 

Tp(z > 1) - y^z^ and Tp{z > 1) - y^z^ (41) 

For non-relativistic fermions, z ^ 1, assume first that the interaction strength y is such that Tp ^ 1 and Tp <C 1. 
Then the interaction terms can be ignored and one recovers the standard equation of state for a free non-relativistic 
gas, a polytrope of 7 = 5/3. For larger values of the interaction strength, the interaction terms become important 
already in the non-relativistic regime, more drastically for the pressure, eq. (|37|) . than for the energy density, eq. H36|l. 
as for the same value of ?/, Tp ~ z^ while Tp ~ z. Hence, the curve for the interacting case will be above the one for the 
non-interacting gas, an effect which increases with increasing interaction strength y. In the region where Tp ^ 1 and 
Tp ^ 1, the interaction terms dominate and the equation of state reduces to p' = p' , which is the stiffest equation of 
state consistent with special relativity. This equation of state was first considered by Zeldovich (see the discussion 
in Note, that this condition is fulfilled already for z ^ y^^^^, which for large values of y is substantially smaller 

than one. Therefore, even for non-relativistic fermions one can reach the causal limit, where p' = p', for an strongly 
interacting gas at p' <C 1. 

For relativistic fermions, z 1, consider first the interaction strength y to be such that Tp 1 and Tp ^ 1. Then 
the free gas terms in eqs. (|38|l and H39I) can be ignored and the equation of state reduces simply to p' = p' . Thus, 
both in the relativistic and in the non-relativistic limit, the equation of state will have the form p' — p' . For small 
values of y, a different regime for the equation of state is reached when Tp <C 1 and Tp <^ 1. In that case the equation 
of state becomes the one of an ultra-relativistic gas p' = p' /3, see eq. (|17|l . This happens for e.g. y = 0.01 and z ~ 10 
in Fig. 21 

With the above arguments we are now in the position to explain the behaviour of the curves in Fig.^] For y ~ 0.01 
and for relativistic fermions, the equation of state is that for an ultra-relativistic free gas, p' = eq. (|17|l which 
differs from p' = p' , the one for an interaction-dominated Fermi gas with large values of the interaction strength 
y ^ 1. The different prefactor explains the slight difference between the lines for y — 0.01 and the other curves with 
larger values of y for large densities p'. There appears a sharp change of the slope of the curves which signals the 
transition from the energy density being dominated by interactions to being dominated by the (free) kinetic terms. 
For even smaller densities p' <C 1, the equations of state are eventually given by the non-relativistic polytrope of 
eq. (|16|l . Deviations from the non- interacting case arise for intermediate densities with increasing interaction strength 
y, since the interaction increases the pressure more rapidly than the energy density. 

2. Solution to the TOV equation for interacting fermions 

We solve the dimensionless TOV, eqs. (|14|) and (|15|) . with the dimensionless pressure and energy density as given 
by eqs. (|34|l and H35() . The equation of state depends now on the interaction strength parameter y. The plot of the 
dimensionless mass M' versus radius R' for different values of y is depicted in Fig. [3 on a double- log scale. 

From Fig. [3 we observe that: 

1. For small interactions strengths y <C 1, the mass-radius curves are very close to each other. The interaction 
terms can be neglected in this case and the equation of state is determined by the free gas term, which results 
in nearly the same mass-radius curves. 
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2. For large interaction strengths, y > 1, the mass-radius curves are shifted towards larger masses and radii. As 
the interaction becomes stronger, the equation of state is getting more and more stiff approaching p' — p' . 

3. The mass-radius curves, in particular the maximum mass and the corresponding radius, seem to follow a common 
trend for strongly interacting fermions, so a scaling behaviour should be expected. 

4. The slope of the mass-radius curve for masses well below the maximum mass changes with increasing interaction 
strength y from the free case behaviour M' x R''^ = const, to one with a rather constant radius for a wide range 
of masses. 

In order to obtain the dimensionful mass and radius for different values of y one should use eqs. H29|l and l|3()|l 
and take the numerical values for Af and R' from the corresponding curve in Fig. |S1 For a neutron star, strong 
interactions between neutrons mediated by vector mesons sets the typical interaction strength to be around y — 10. 
Then, the maximum stable mass emerges to be 3.77Mq instead of O.TIMq as was found earlier for non-interacting 
neutrons. Note, that attractive forces will reduce this value somewhat closer to the maximum masses considered 
presently for realistic equations of state (see e.g. For an interacting neutralino star, weak interactions result in 

low values for the interaction strength, y ~ 1/3, so that the maximum mass would stay unchanged. However, if we 
assume strong interaction values for y as motivated by QCD of around y — 10'^, the maximum stable mass is boosted 
to 2.7 X 1O~^M0 compared to 6.3 x 1Q~^Mq for the non-interacting case. 

Next we discuss the behaviour of the mass-radius relation for large values of R' considering strongly interacting 
fermions, i.e. y 3> 1. The question is, is there any relation similar to the one for the non-interacting case, eq. H27|l ? 
Note first, that the tail of the M' versus R' curve corresponds to small central densities of the star p'(0), so that the 
dimensionless Fermi momentum z = kf/nif is quite small, i.e. z <C 1. In the non-relativistic regime, it is likely that 
the equation of state can be described by a polytrope, even for interacting fermions. 

Indeed, at z <C 1, the equation of state reduces to eqs. H36|l and H37I) and one can ignore the interaction term 
proportional to as compared to the term proportional to z^. Thus one finds that 



(42) 
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Now, using this equation to eliminate z from eq. 1)3 7|) . we get 



(43) 

For large values of y 3> 1, the interaction term becomes important for the pressure above a certain Fermi momentum, 
when rp(z <C 1) '^-^ 1 or z ~ l/?/^- Then, the pressure grows with the square of the energy density and the corresponding 
polytropic coefficient is 7 = 2. Using the Lane-Emden solution, (|28|l . one finds that the radius does not depend on 
the central density, while the mass increases linearly with the central density, hence i?' = const, for a large range of 
masses. Examining Fig. ^ we see that this behaviour is present for ?/ 3> 1 close to the point where the causal limit, 
p' — p' , is reached. The upper limit to the energy density is therefore given by the point where the interaction starts 
to affect the energy density in addition to the pressure, i.e. Tp{z <^ I) ^ 1 or p ^ ^ y~^. The critical lower value 
is given as discussed above by the condition Tp{z <C 1) ~ 1, i.e. z ^ or p' cx y~^, where the interaction starts to 
affect the pressure. Therefore, radii are constant for the mass range extending from M' oc p' oc to M' cx p' cx 
which can reach, for the case y = 10'^, up to about twelve orders of magnitude in compact star masses! 



3. Scaling behaviour for interacting fermions 

As mentioned above, a scaling behaviour is expected for large values of y. Consider eq. H36() and eq. (|37|l for non- 
relativistic fermions and for y 3> 1. Then the dimensionless pressure p' is dominated by the interaction term while 
the dimensionless energy density p' , (|36|l . contains both the kinetic and interaction terms. One can see that in this 
case the equation of state will turn from the non-relativistic polytrope, (|16|) . to the causal limit of the form p' — p' 
when the fermions are still non-relativistic, i.e. z <^ I. Indeed, one can write the dimensionless energy density as: 



z3 1 



z 
3^ 



3 



p' (44) 



Hence, the pressure p' is larger than the kinetic term of the energy density for z^ ^ y^z^ or z^ ~ 1/?/^- In this regime, 
both p' and p' are of the order of ~ Hence, the pressure and energy density can be rescaled in a dimensionless 

form with the factor m'^/y'^. The corresponding Landau mass and Landau radius will be modified accordingly to 

Mr'=M^/m}-y and RT'^Mp/mj-y , (45) 

respectively, so that the maximum mass and the corresponding radius increase linearly with the interaction strength 
y. We note in passing that the importance of the interaction on the global properties of compact objects has been 
also noted for the case of boson stars with interacting scalar fields in [53 . 

The numerical results for the scaling behaviour of the maximum mass as a function of the interaction strength y 
are plotted in Fig. on a double-log scale. For small values of y, the maximum mass M^^^ does not change and 
basically remains constant. For y 1, the maximum mass M'^^^^ starts increasing as a function of y, approaching a 
power-law rise for y 3> 1. 

For small values of y, the interactions do not affect significantly the equation of state and there will be almost no 
change in the maximum mass. If the interaction terms in the equation of state become comparable to the ones of the 
free gas, i.e. if y ~ 1, the maximum mass begins to increase. Finally, for y 1, the pressure and the energy density 
for the maximum mass configuration are dominated by the interaction terms, which can be scaled out by the factor 
1/y^ so that the maximum mass will increase linearly with y as explained above. For the part of the graph in Fig.|H| 
where y 3> 1, one can perform the following general fit: 

M;„a. = ci + si • , (46) 

as log M^jj^ is linear in log y for y 3> 1 . The constant ci is fixed by the numerical result found for the non- interacting 
case, hence c = 0.384. The parameters si and 71 are fitted to the curve of Fig.|Hlto be si — 0.165 and 71 = 0.999 « 1. 
The maximum mass Mmax of a compact star can then be calculated from the relation 

Mmax = (0.384 + 0.165 • y) ■ (^-^^\ ■ 1.632Mq (47) 

V "1/ y 



as an approximation to our numerical solution of the full TOV equations. 
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FIG. 6: The dimensionless maximum mass M^^^ for interacting fermions versus the interaction strength j/ on a double-log 
scale. For small values of y, the maximum mass does not change, while it increases with a power law for strong interactions 

(y > 1)- 



DC 



1000- 



100 



10- 



Weak 
Interaction 



Strong 
Interaction / 



1 

0.01 0.1 



m, I c 

100 1000 



1 10 

Interaction Variable y 

FIG. 7: The minimum dimensionless radius Ji^jn for the maximum mass configurations as a function of the interaction strength 
?/ on a double-log scale. 



16 




We plot the minimum radius i?J„j„ versus y on a double-log scale in order to extract a scaling behaviour for the 
radius in Fig. [7| For small values of y, the radius stays constant. For y > 1, the radius increases with the 

interaction strength following again a power-law behaviour, so that we take the following fitting expression for R'^i^: 

Kmn - C2 + S2 • y^^ (48) 

where C2, S2 and 72 are constants. The constant C2 is fixed by the non- interacting limit as discussed before, i.e. 
by i?^m for a free Fermi gas, for which numerically C2 — 3.367. In the y S> 1 limit, one obtains S2 = 0.797 and 
72 — 0.9942 ~ 1. Hence, for large interactions strengths y ^ 1 the radius can be approximated as 

i?^i„ = (3.367 -I- 0.797 -y) • ( ) •2.410 km (49) 

V "1/ / 

According to the scaling arguments as derived above, the minimum radius (as well as the maximum mass) should 
increase linearly with y for large interaction strengths which is indeed being found numerically and observed in Fig. El 
for the maximum mass and in Fig. [7| for the corresponding radius. For small values of y, scaling arguments predict a 
nearly constant minimum radius as a function of y as also clearly seen in Fig. and in Fig. [3 

However, one should expect a difference in the radial density profile for the weak and the strong interaction cases. 
When interactions are weak, the equation of state is governed by the ideal Fermi gas terms in the core of the compact 
star. For increasing interaction strength, the pressure increases more rapidly due to the interaction terms for low 
energy densities. As a result, at the surface of the compact star, the energy density decreases more rapidly and the 
additional pressure from the interactions becomes important. A halo is created as more matter can be supported 
against the gravitational pull. The total mass of the star remains almost the same as it is mainly given by the dense 
core and not affected by the dilute halo. Interestingly, also the corresponding radius is not modified substantially by 
this effect. This behaviour is illustrated in Fig. |H1 which shows the density profile for various star configurations with 
different values of the interaction strength y. The dimensionless density p'{r') is plotted versus the dimensionless 
radius of the star r' for several configurations. Note, that the curves are drawn for the maximum mass configurations 
and normalised to their central density and the total radius. Fig. |S1 shows the presence of a dense core and a low 
density halo at the surface of the star. The rate of the decrease of p'(r') with the radius r' is large for small values of 
y and decreases gradually for higher values of y generating a denser halo of the compact star for strongly interacting 
fermions compared to the weakly or non-interacting ones. 
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V. SUMMARY AND DISCUSSIONS 

We construct general equations of state for fermions of arbitrary mass rn/ to be used as input for solving the TOV 
equations for corresponding compact stars. Besides a free gas of fermions, wc consider also the case of interacting 
fermions by adding interaction terms ~ to the energy density and pressure. We discuss various ways of rescaling 
the TOV equation and find the corresponding scaling solutions for compact star configurations for arbitrary fermion 
mass and interaction strength. The scaling solutions are tested by detailed numerical calculations. 

In particular, we have demonstrated that Landau's argument and the corresponding expressions for the maximum 
mass and the corresponding (minimum) radius hold also in the case of general relativity. For a compact star made 
of fermions, there exists an upper limit for the mass which is of the order of the Landau mass Mp/rnj-, where Mp 
is the Planck mass. The corresponding radius scales as the Landau radius Mp/rnj:. By directly solving the TOV 
equations, we have found the mass-radius relation in dimensionless form, i.e. in units of the Landau mass and Landau 
radius, supporting our analytic finding. The numerical values of the maximum mass and the corresponding radius 
are M^ax = 0.384Mp/m'j and Rmin = 3.367 Mp/mj, respectively. For large radii the mass and the radius are related 
by the well-known relation M ■ B? = const. The results are in accordance with those known for a free neutron gas 
but our scaling relations can be used for fermions of any mass. For example, fermions with a mass of 100 GeV can 
form compact objects up to a maximum mass of about 1O~^M0 and radii down to 1 meter as deduced from the same 
mass-radius relation as for a free gas of neutrons, when properly rescaled by the Landau mass and Landau radius. 
The same can be done for the mass-radius relation of neutrino stars. For m/ ~ 1 eV the corresponding maximum 
mass and its radius are, however, of the order of Mmax ^ lO^^M© and Rmin ^ 10^^ km. Interestingly, a hypothetical 
fermion with a mass of about 10^^ eV can have a maximum mass and a corresponding radius which matches the total 
mass and the horizon of the present universe. 

Effects from interactions between the fermions are taken into account by adding terms proportional to the fermion 
density squared to the expressions for the pressure and energy density. This can be motivated by e.g. such effective 
models of strong interactions as quantum hadrodynamics. The resulting equation of state depends on a new variable 
y = mf/rnj which controls the strength of the interaction. The mass mj fixes the range and strength of the 
interactions. The values of y can be as small as 10~^^ for neutrinos with a mass of ~ 1 eV and as high as y = 10'^ for 
strongly interacting neutralinos with a mass of ~ 100 GeV. 

Using scaling arguments we arrive at the dimensionless equation of state and corresponding scaling solutions for the 
maximum mass and its radius for interacting fermions. We show that the maximum mass and the radius are rather 
constant for small interaction strengths, y <C 1, but are entirely determined by the interaction terms for y » 1: Both, 
the maximum mass and its radius increase linearly with the interaction strength for y ^ 1. The scaling is supported 
by numerical calculations, where we find that M^ax ~ 0.165y- Mp/mj and Rmax ~ Q.797y- Mp/mj for y ^ 1. Besides 
these general scaling features, there is a small change in the density profile for large interaction strengths, which can 
be attributed to the formation of an enhanced halo density in the outer region of the star. For large interaction 
strengths the mass-radius relation changes from the standard M ■ R"^ = const behaviour to the one with a constant 
radius for a wide range of compact star masses. This is because the pressure increases as the square of the energy 
density and not as the power 5/3 for a free non-relativistic Fermi gas. The mass range for constant radii changes 
with the interaction strength in the range between y~^ and y'^. Specifically, for values of the interaction strength of 
y = 10^ and a fermion mass of 100 GeV, the maximum mass increases from the value for the non-interacting case, 
Mmax ~ 1O~^M0 to Mmax ~ 1O~^M0, which is comparable to the one for ordinary neutron stars. The compact 
star mass range, where the radius stays constant, extends from about the maximum mass of about 1O~^M0 down to 
1O~"'^^A^0 with a typical radius as given by the Landau radius, i.e. about 100 meter. 

There are two important issues which require a special study which is beyond the scope of the present paper. 

First, the crucial assumptions for all these investigations is that the fermions are stable on the time scale comparable 
with the lifetime of the universe, i.e. r > H^^ w 14 Gyr, where Ho is the present value of the Hubble constant. In 
other words, it is assumed that the fermions constituting the compact star are conserved, i.e. there is no annihilation 
into other kinds of matter. Naive estimates for the lifetime using r (n • with the number density n 

and the cross section a ^ mj/mj, result in the constraint nij < (Hq ■ m.jy/^. For weak interactions, one arrives at 
fermion masses w/ < 1 keV, for gravitational interactions at the Planck scale mi = Mp, however, the fermion mass 
must be only lower than m/ < 10^ TeV (in all cases y <^1). 

The second important question is, when and how the compact objects made of exotic fermions could be formed? 
One may speculate that this could happen at very early stages in the history of the universe, right after the inflation 
stage. These early formed objects could serve as seeds for clumping ordinary matter at later stages of the expansion, 
after the radiation decoupling. Therefore, one may speculate about hybrid objects where exotic fermion stars are 
surrounded by the halos of ordinary matter. 

Our final remark concerns the possible observations of compact objects made of dark matter particles. In fact, 
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there exist limits derived from the observation of micro-lensing events. The MACHO collaboration has excluded the 
mass range of (10~^ — 30)Mq for compact objects forming the bulk of the Galactic dark matter. However, compact 
stars of these mass ranges are not ruled out if they do not contribute more than 4 • IQ^^Mq to the Galactic halo (see 
[56| and references therein). 
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